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Analyses and experiments demonstrate the potential benefits of optimizing piston and displacer motion in a 
free-piston Stirling Engine. Isothermal analysis shows the theoretical limits of power density improvement due 
to ideal motion in ideal Stirling engines. More realistic models based on nodal analysis show that ideal piston 
and displacer waveforms are not optimal, often producing less power than engines that use sinusoidal piston 
and displacer motion. Constrained optimization using nodal analysis predicts that Stirling engine power density 
can be increased by as much as 58% using optimized higher harmonic piston and displacer motion. An 
experiment is conducted in which an engine designed for sinusoidal motion is forced to operate with both second 
and third harmonics, resulting in a piston power increase of as much as 14%. Analytical predictions are 
compared to experimental data and show close agreement with indirect thermodynamic power calculations, but 
poor agreement with direct electrical power measurements. 


I. INTRODUCTION 


TIRLING engines can achieve high thermal efficiency (net work out/heat in) but have low power density (power 

per unit volume) when compared to many open-cycle engines'?*. One potential method of increasing power 
density is to enforce piston and/or displacer motion that more closely approximate those of the ideal Stirling cycle. 
Achieving the ideal cycle requires that the piston and displacer dwell and abruptly change direction throughout the 
cycle, which is difficult to achieve in both kinematic and free-piston configurations. Designers typically sacrifice ideal 
motion for practically achievable motion, most often converging on a mechanical linkage or electrical control scheme 
that imposes sinusoidal or nearly sinusoidal motion*®. While there are certainly benefits to choosing sinusoidal 
motion, it is not the optimal choice since sinusoidal motion reduces power density. There are some patents based on 
achieving piston motion that more closely approximates ideal motion in kinematic Stirling engines”*®, but the author 
has not found any examples of a functional kinematic Stirling which uses non-sinusoidal motion for the purpose of 
improving performance and no examples of attempts to achieve performance benefits through non-sinusoidal motion 
on free-piston Stirling engines. 

Stirling cycle applications include cryo-coolers*"®, natural gas co-generation units , solar-dynamic power 
conversion!*!>'6, and nuclear dynamic power conversion!”!®:!9:22!_ They are typically used in applications which have 
high fuel costs or in systems that require closed-cycle operation. High efficiency and closed-cycle operation are both 
requirements of space power systems, making free-piston Stirling engines candidates for these applications. They are 
a key technology in NASA’s Radioisotope Power System program because their high efficiency potentially enables 
NASA to increase the number of missions it can fly over the coming decades using the limited supply of plutonium- 
238°9*12. They also trade favorably in low to moderate power level fission power applications because their high 
efficiency requires less heat input from the reactor and reduces heat rejection requirements, reducing the mass of the 
reactor shield and the radiators'”'®'°.7°. Stirling engines have been considered for use in several terrestrial applications 
including automotive engines, concentrated solar power plants, and residential co-generation systems, especially when 
rising fossil fuel costs indicate that their high efficiency can make them economical. However, the low power density 
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and high cost of Stirling engines typically make them unpractical when fuel costs are low and internal combustion or 
other open-cycle engines are viable. 

Figure 1 shows a schematic of the ideal piston and displacer position for a Stirling engine in the beta configuration. 
Process 1-2 shows constant temperature compression, 
in which the displacer remains still, minimizing the 
expansion space volume. The piston compresses the 
gas as heat is removed through the cooler. Process 2- 
3 shows constant volume heat addition, in which the | 
displacer moves gas from the cold side to the hot side, © 4 
through the regenerator, while the piston remains still. | 
Process 3-4 is a constant temperature expansion Time 
process, in which the expansion of the gas moves the Figure. 1. Piston and Displacer Waveforms for the Ideal Stirling 
piston, and the displacer moves along with it to Cycle 
minimize volume of the compression space. The work 
done during the expansion process minus the work done during the compression process is the usable energy produced 
by the Stirling cycle. In free-piston engines this power is typically extracted through the linear alternator. Process 4- 
1 is a constant volume heat removal process in which the displacer moves gas from the hot to the cold side, through 
the regenerator, while the piston remains still. For an ideal Stirling engine, as the regenerator effectiveness approaches 
unity the external heat addition and rejection requirement during the constant volume processes approaches zero, and 
the Stirling cycle efficiency approaches the Carnot efficiency. The power density of an ideal Stirling engine with zero 
dead volume is independent of regenerator effectiveness and is proportional to the natural log of the volume ratio. 

Ideal piston and displacer motion requires the piston and displacer to dwell during periods within a cycle and abruptly 
change direction during others, requiring large forces. The commonly used alternative to ideal motion is sinusoidal 
motion, which can be analyzed under isothermal 
assumptions using the method proposed by Gustav Specific Power Ratio for t= 2.5 
Schmidt‘. Isothermal engines with perfect 
regeneration operating under the Schmidt cycle —— Esta Leis af Piston Motion 
achieve the same efficiency as the ideal Stirling 
engine, which is equal to the Carnot efficiency. The 
sinusoidal piston motion associated with the Schmidt 
cycle, however, does reduce power density. 

Isothermal analysis of the ideal and Schmidt cycles = = 
shows the extent to which theoretical power density —~ 
can be improved through the use of ideal motion. 
Figure 2 shows the ratio of ideal to Schmidt cycle 
power assuming equal limits of piston motion or 
equal working space amplitudes. The complete Figure 2. Comparison of power density ratio for ideal and Schmidt 
analysis is shown in Briggs? with major points cycles assuming equal limits of piston motion or equal working 


Rat 


highlighted below. space amplitude 
-_ r+Kq(t Tq-1) - r +Kq(Tq—-1) 
Wnet,ideal ~ mRT cs [z In (ae) In (eC ) 
= Vinax = Tes = Vp ee TesCy 
Vinita Tes : Vinita . Vp 


V, V, V, V, V, 
fs ( esti Veg UR Way eset) 
Tos Tc Tepper Tu Tes 


1-v1-b2 , 
Wscnmiat = 2™MRT cs (=) (=) (t —1)sina(2) 


Vt? + 2Ktcos(a@) + K2 
Kr 


_ Ves,sw _ 2TgsCp 
ees =o 


Vessw Vessw 


Cz = 


(“eee Vesct Vo Vr Vir Vescl + “sexe 
2Tcs_ Tes = ToS OTrpgr TH =Trs Ts 

Where m is the working gas mass, R is the specific gas constant, T is the temperature of a given gas space, and Vis 
the volume of a given gas space. The subscript max refers to the maximum cycle value, min refers to the minimum 
cycle value, ES refers to the expansion space value, CS refers to the compression space value, D refers to the dead 
space value, C refers to the cooler space value, CL refers to clearance volumes, r refers to regenerator values, H refers 
to heater values, and SW refers to swept volumes 


II. ANALYSIS 


A. Nodal Analysis - Sage Stirling analysis software 


Isothermal analysis assumes that heat transfer rates in the expansion and compression spaces are so high that the 
temperature of the gas within those volumes can be considered isothermal. Real engines lack sufficient surface area 
in the compression and expansion spaces to accomplish this, so high surface area heaters/acceptors and 
coolers/rejectors are added along with their associated dead volume. 

Improving model accuracy beyond what is achievable with isothermal or adiabatic analysis requires nodal analysis 
and a specific engine design, including heat exchanger geometry, regenerator specifications, cylinder geometry, seal 
geometry, etc. One such nodal analysis tool is a commercially available program called Sage*®. Sage is commonly 
used for Stirling engine optimization in the design phase. The Sage is a one-dimensional, cyclic steady-state model 
that couples the equations of motion of the piston and displacer with the Navier-Stokes equations and energy equations. 
Sage can also be used to calculate the reduction in available energy, which is ignored by isothermal and adiabatic 
analysis. 

Sage is used to model a 1-kW Stirling engine which is representative of engines tested at GRC. This model assumes 
an isothermal boundary condition on a solid surface node within the engine. it then calculates temperature gradients in 
the rest of the solid and the gas based on solid conduction and calculated convection coefficients. The axial temperature 
distribution along the base of the finned exchangers (acceptor and rejector) are set as an input. These temperature 
inputs typically come from heat transfer analysis done outside of Sage. Sage then iteratively solves for gas temperatures 
and pressure drops by guessing and checking fin temperature profiles and gas velocities and the resulting displacer 
motion (piston motion is typically a user input). 

The predicted power output for isothermal analysis is 36% larger than the Sage predicted power output. Isothermal 
analysis predicts an engine efficiency equal to the Carnot efficiency, which in this case is 0.550 (Tn = 779 °C and T; = 
350 °C) which is 80% higher than the efficiency of 0.307 predicted by Sage, suggesting that isothermal analysis does 
a relatively poor job of predicting the heat transfer requirements in real engines. 

Sage requires all non-sinusoidal time dependent inputs to be entered as sums of harmonic sinusoids. As higher 
harmonics are added, the solver time resolution must be increased to resolve the higher frequencies, increasing 
computational time. Ideal piston and displacer motion are approximated using 7-term truncated Fourier series (R? 
values for all cases run are above 0.995). 

One inherent problem in modeling ideal piston and displacer motion in Sage or any other high fidelity Stirling model 
is that there are in infinite combination of piston and displacer motions which are considered ideal. Piston motion can 
be symmetric or asymmetric, with varying dwell times and many combinations of dwells or piston/displacer constraint 
methods are still considered ideal. Since isothermal analysis assumes that heat transfer takes place instantaneously and 
does not take into account pressure drop or other loss mechanisms that are functions of piston and displacer velocity, 
predicted power output is equal for all ideal waveforms and the efficiencies each of these waveforms are equal to the 
Carnot efficiency, regardless of dwell times or piston/displacer velocities. However, higher fidelity tools such as Sage 
consider the effect of insufficient dwell times on gas temperatures and higher mass flow rates on pressure drop through 
the regenerator and heat exchangers, so they predict that different versions of the ideal waveform have different power 


output and efficiency. 


B. Nodal Analysis — Comparison of Sinusoidal and Ideal Waveforms 


In this section nodal analysis is used to compare performance of ideal waveforms with sinusoidal waveforms. The 
ideal motion examined in this section requires external forcing on the piston and displacer. In reality, most free-piston 
Stirling engines have a free displacer attached to a linear spring which resonates near the convertor operating frequency. 
These engines have no mechanism to impose an ideal waveform on the displacer. However, kinematic Stirling engines 
and some unproven free-piston designs do have mechanisms to impose non-sinusoidal displacer waveforms. The 
analysis in this section is intended to show the effect of ideal piston and displacer motion on engines designed with 
such capability. Engines operating with non-sinusoidal piston and displacer motion are compared to an engine 
operating with a sinusoidal piston and a free displacer. 

Ideal piston and displacer motion result in maximum power density and efficiency in ideal Stirling engines. In real 
engines, increased piston and displacer velocities inherent to ideal waveforms can adversely affect engine performance 
by increasing pressure drop through heat exchangers and/or adversely affecting heat transfer. The following analysis 
shows that these parasitic effects overwhelm the increase in thermodynamic power provided by ideal waveforms when 
applied to this engine geometry. 

This analysis constrains the maximum and minimum working volume, minimum expansion space volume, and 
minimum compression space volume to the same values as sinusoidal motion. These constraints limit the inward travel 
of the piston and displacer, leaving a stagnant gas volume in the compression space and limiting the working space 
stroke. Figure 3 shows a comparison of ideal and 
sinusoidal motion under these constraints. Isothermal 
analysis predicts that under these constraints the 
power density of an engine operating with ideal 
piston and displacer motion is 50% higher than the 
same engine operating with sinusoidal piston and 
displacer motion, with no change in efficiency Ref 
25. Sage-based nodal analysis of ideal motion under 
these constraints is more pessimistic, predicting 
lower power density and efficiency than predicted for 
sinusoidal motion in most cases. 

Four versions of the ideal Stirling piston/displacer Figure 3. Comparison of ideal and sinusoidal motion assuming 
waveform are selected for study, three symmetric equal minimum spacing between piston and displacer. 
piston waveforms with varying piston dwell times and one asymmetric waveform in which the piston dwell times are 
selected by minimizing the piston and displacer RMS velocity in an effort to reduce pressure drop losses through the 
regenerator. The results of nodal analysis applied to these waveforms appears in Table I. 


Regardless of dwell time, symmetric ideal waveforms TABLEL 

operating under the constraints of Case 1 produce less NODAL ANALYSIS RESULTS FOR IDEAL WAVEFORMS 
power and are less efficient than an engine operating with Net | Piston | Displacer | Efficiency 
sinusoidal piston motion and free displacer motion. This is POWER 20m et) « AONE (%) 
due to the increased viscous dissipation in the regenerator - - (Ww) Ow (Ww) 

; ; . Sinusoidal 1170 1170 0 30.7 
and heat exchangers caused by increasing piston and/or Ssumne inc 
displacer velocity. The medium dwell case is the bestof the | 4.44 4 50 1680 1630 12 
three symmetric cases examined, but it offers no Symmetric 
improvement over sinusoidal motion. Removing the | dwell=7/2 820 |) 1580 760 16.5 
symmetry constraint on ideal waveforms allows the piston | Symmetric 
and displacer to achieve ideal motion with lower velocities, | dwell = 37/4 630 1320 690 12.7 
allowing for improved engine performance. This isa coarse Asymmetric i150 i580 Aso 218 
analysis of the dependence of engine performance on piston | Min Velocity 


dwell time, and it is likely that an optimum performance 

point exists at dwell time between the points chosen for analysis above. It is possible that some version of the ideal 
waveform could show an increase in power density over sinusoidal motion, but the results of this analysis suggest that 
replacing sinusoidal waveforms with ideal waveforms most likely does not produce a combination of power and 


efficiency that sufficiently improves on sinusoidal motion to justify implementation. 


C. Optimized Waveforms 


It may seem discouraging that under the chosen constraints there is no benefit to achieving ideal motion for this 
engine design. However, seeing that ideal motion results in increased piston power in all cases suggests that there may 
exist some optimal motion, different from ideal motion, that takes advantage of the increased piston power without 
incurring large penalties due to increased viscous losses through the regenerator and heat exchangers. Determining 
this optimum motion requires the use of non-linear constrained optimization. This is achieved using the fmincon 
function in MATLAB, using net power, as predicted by 
Sage, as the objective function. Various non-linear 
constraints can be used depending on how one wants to 
match sinusoidal and optimal waveforms. In all cases the 
out-limits of piston and displacer motion are set equal to the 
values generated in the sinusoidal case, guaranteeing that 
the engine does not have to grow to accommodate ideal 
motion. The following cases are studied in which different 
characteristics of the sinusoidal waveform are used as Angie radians} 
optimization constraints: 


1. In-limit of piston equal, minimum piston/displacer 4.000 
gap equal 

2. In-limit of piston equal, no limit on displacer inward = 
travel é ; : 

3. No constraint on the in-limit of piston or displacer oe, eA eee 
motion 
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The optimization uses a seven term Fourier series to —Sinuzoith [PetonFDdagam —— -4.000 
define the piston and displacer waveforms, giving 27 
independent variables (seven piston amplitudes, six piston 
phase angles with the phase angle of the fundamental 
frequency is fixed at zero, seven displacer amplitudes, 
seven displacer phase angles). Using single term sinusoidal 
motion as starting point, optimization requires hundreds of 
iterations and thousands of function evaluations. Each 
function evaluation requires a converged Sage solution, 
making the optimization a computationally intensive cei eens 300 
process. The optimization results are summarized in Table 
Il. 

In each case the optimized waveform produces more 
power at reduced efficiency than the sinusoidal case, 
opening a viable trade space between power and efficiency for engine designers. Figure 4 compares sinusoidal and 
optimized piston/displacer motion, piston force displacement diagram, and displacer force displacement diagram for 
Case 2. The complete dataset can be found in Ref 25. 

This optimization uses only power output in the objective function. An alternative would be to optimize on a 
weighted combination of power and efficiency in which the weights are chosen according to the specific application., 
however a full mapping of the power/efficiency trade space 
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Figure. 4. Comparison of Sage results for sinusoidal and 
optimized waveforms under Case 2 constraints. 


i : TABLE II 
is, as yet, unstudied. NODAL ANALYSIS RESULTS FOR OPTIMIZED WAVEFORMS 
Net Piston | Displacer | Efficiency 
III. EXPERIMENTAL TESTING Power | Pow | Paws (%) 
The numerical analysis above suggests that Stirling (W) (W) (W) 
engine power density and efficiency can be traded against | Sinusoidal | 1170 | 1170 | 0 | 30.7 | 
each other and specified to the designer’s needs. This Case 1 1430 1580 150 27.9 
section describes a test performed in the Stirling Research [Case2_ | 1630 | 1730 | 100 | 66 
Case 3 1850 1853 3 25.3 


Lab (SRL) at the NASA Glenn Research Center intended to verify that non-sinusoidal piston motion is achievable 
and to verify the predictions of nodal analysis. 


The engines chosen for testing are a pair of 1-kW engines previously used to test a sodium-potassium (Nak) heat 
exchanger’’ and to test the feasibility of operating two engines with a common expansion space*® for the Fission 
Surface Power project. The specific cases studied above could not be replicated with these engines because they 
have no method of altering the displacer waveform. Therefore only the piston waveform is altered and the displacer 
freely responds to the resulting changes in the pressure wave. In practice the stiff linear spring attached to the 
displacer results in nearly sinusoidal displacer motion in all cases. 


Since the piston is not designed to resonate at higher harmonics, higher harmonic motion had to be imposed by the 
alternator, increasing the spring force required from the alternator and the alternator current. The current required to 
run non-sinusoidal piston waveforms at the nominal operating amplitude of 10 mm saturates the alternator before the 
waveforms can be substantially altered. However, measurable changes in piston motion and engine performance can 
be achieved when operating at reduced amplitude. 6 mm is chosen as the baseline for the experiment because it is 
the minimum recommended amplitude required to charge the internal gas bearings. Alternatively, if an engine were 
available that was designed for multi-harmonic resonance the alternator force would be reduced and greater 
amplitudes could be achieved and alternator losses could be minimized. 


The out-limit of piston motion for all tests is set equal to the value measured during sinusoidal testing. The first 
harmonic of the piston amplitude is then increased and additional harmonics are added until the out-limit of piston 
motion returns to the value measured during 6 mm baseline testing. The non-sinusoidal piston waveforms are 
allowed to exceed the in-limit of piston motion, consistent with constraint Case 3. 


Non-sinusoidal piston waveforms can be obtained using a programmable AC source to impose an arbitrary voltage 
waveform on the alternator terminals. However, the current requirements of this method is prohibitive. Instead, 
three AC sources are connected in series, each dedicated to a single harmonic and paralleled to a capacitor bank. 
This method reduces the current requirement of the AC sources, but also limits the number of harmonics that can be 
tested due to manufacturer limits on the number of AC sources that can be put in series. 

Alternator inductance causes the voltage measured at the alternator terminals to be larger than the EMF and out of 
phase with the alternator current (even when the engines are operated at resonance). In single harmonic engines, a 
tuning capacitor can be used to compensate for the alternator inductance so that the controller voltage is reduced to 
the level of the EMF and the current and “post-capacitor” voltage are once again in phase. However, when additional 
harmonics are introduced, the tuning capacitor cannot cancel the alternator inductance of the higher harmonics so the 
voltage required of the control system is much higher than the EMF seen by the piston. There may be an alternative 
solution to this problem involving changes to the control system electronics, but this experiment used the “brute 
force” approach, in which the AC sources used as the control system supply relatively large voltages in order to 
impose relatively small piston EMF. For example, a peak-to-peak third harmonic voltage of 420 V results in an 
estimated peak-peak EMF of 30 V (estimated EMF based on piston velocity measurements times the motor 
constant). 

The methods used in this experiment are not ideal for demonstrating system level benefits of non-sinusoidal 
waveforms on Stirling engines. However, the intent of this experiment is to achieve a proof-of-concept on existing 
and available hardware and provide an experimental database for comparison to model predictions. Sophistication 
and optimization of the overall engine and control system is left as future work. 

Heat input is provided using cartridge heaters capable of delivering 4 kW of power to each engine. Voltage to 
each of these heaters is regulated using a Variac. Each engine uses a separate Variac so that the average hot-end 
temperatures of the engines can be made equal despite heater cartridge variability and differences in thermal 
resistance between the cartridges and the engines. Heater power is measured and is used to calculate gross engine 
efficiency. The cold end of the engine is water cooled using a recirculating water chiller. 

Piston and displacer motion are measured using Fast Linear Displacement Transducers (FLDT). Alternator current 
is measured using a Pearson coil on each engine. The FLDT signals are read into an oscilloscope to capture the 
waveforms in real time. A high precision power meter, which is also capable of recording waveforms, is used to 
measure alternator current and pre-capacitor voltage. The post-capacitor voltage and alternator current are read by a 
separate oscilloscope. Both oscilloscopes and the power meter are synchronized using the piston FLDT signal to 
guarantee proper phasing during post-processing of the data. 


Hot-end temperature on the engine is measured using eight type K thermocouples. Cooling water temperature on 
each engine is measured using eight platinum RTDs (four each on inlet and outlet). Alternator housing temperature 
and bounce space gas temperature are also measured using four thermocouples on each engine. 

Steady-state measurements are taken at a rate of 1 Hz. Higher speed data, on the order of kHz are taken using the 
oscilloscopes and power meter described above. 


A. Methodology 


Each test begins by establishing a steady-state baseline at 6 mm amplitude, 550 °C hot-end temperature and 50 °C 
Cold-end temperature. The baseline point uses only the first harmonic power supply, producing sinusoidal piston 
motion, which is the typical operating mode of free-piston Stirling engines. This point is used to establish both the 
baseline power measurement and the out-limit of piston motion which is matched by all non-sinusoidal waveforms. 
At the conclusion of baseline testing the AC Bus voltage of the second or third harmonic is increased by 45 V at a 
phase angle of 180 ° relative to the 50 Hz AC Bus voltage. This results in a decrease in the out-limit of piston motion 
as the bottom of the waveform flattens. The out-limit of piston motion is then returned to the baseline value by 
increasing the voltage of the 50 Hz AC Bus. Heater power is then increased to maintain a hot-end temperature of 
550 °C. This process is then repeated until some limit is reached. While the AC Bus increment of 45 V results in a 
relatively coarse test matrix, it is chosen because this increment results in a measurable change in the out-limit of 
piston motion for both the second and third harmonics. Smaller increments could be used if the piston position 
measurement was made more precise and repeatable from cycle to cycle. Throughout testing the baseline point is 
revisited to establish repeatability and to verify that the lower limit of piston motion has not changed. Numerical 
optimization, using the Sage nodal analysis tool, identifies phase angles of -1/2 for the second harmonic and zero for 
the third harmonic as optimizing power output. These are the only phase angles tested in this test experiment. 
Exploration of phase angle sensitivity and influence on other parameters such as efficiency is left to future work. 

Although the control system is capable of providing 850 Vrms to the alternator, engine or alternator limits are 
reached well before these voltages are achieved. Testing beyond 160 Vrms on the second harmonic results in 
unbalanced current between engines, erratic power measurement, and eventually an overstroke protection circuit trip, 
resulting in engine stall. This could be an indication of alternator saturation, but it is well below the engine 
manufacturer’s prescribed limits on the alternator voltage and current. This instability is repeatable and limits the 
second harmonic test matrix to 160 Vrms. Testing beyond 295 Vrms on the third harmonic results in an audible and 
abnormal noise coming from the engine. This noise is not consistent with any type of impact, whether piston- 
displacer or piston-endstop, and could be benign. However, testing beyond 295 Vrms on the third harmonic is 
avoided in an effort to preserve the engines. When operating with both the second and third harmonics 
simultaneously, a similar noise is heard when the third harmonic voltage is increased beyond 105 Vrms. These 
engine/alternator/control constraints established the operational limits for non-sinusoidal testing. 

This experiment aims to compare the thermodynamic power output of free-piston engines operating with several 
different piston waveforms. The engines could not be modified for dynamic pressure measurement within budget 
and time constraints, so it is not possible to produce true P-V diagrams based on direct measurements. However, 
there are several methods of calculating thermodynamic power output based on measured parameters such as 
alternator power, alternator voltage, alternator current, alternator resistance, and piston position. The most direct 
power measurement is taken by time averaging of the instantaneous alternator voltage (pre-capacitor) multiplied by 
the instantaneous current at the alternator terminals. This measurement is made using the internal algorithms of the 
power meter as well as through post-processing of the voltage and current waveforms. Adding the resistive losses of 
the alternator to this measurement gives an approximation of the thermodynamic power of the engine. 

This method is useful for determining average power output of the engines, and is used as the figure of merit in the 
results section. However, due to the alternator inductance it does not give good information on the instantaneous 
piston power or forces, making it less useful for the purposes of model comparison. Instantaneous power can be 
calculated by multiplying the measured current by the EMF calculated by multiplying the piston velocity by the 
motor constant. The time average of this quantity gives an alternate measure of thermodynamic power output of the 
engine. 


B. Results 


The operating conditions for the baseline case are 550 
°C hot-end temperature, 50 °C cold-end temperature, 


and 6 mm piston amplitude. The out-limit of piston 
motion under these conditions is 6.06 mm. The 
measured out-limit of piston motion is kept constant 
throughout higher harmonic testing. Table II 
summarizes the various power measurements and 


calculations described in the Methodology section as 
well as efficiency. The gross efficiency measured 


during testing is calculated by dividing the piston 


power (alternator power plus resistive losses) by the 


gross heat input (electric power supplied to the 


cartridge heaters). The gross heat input includes heat 


that is lost through the engine insulation package, so it 


underpredicts thermodynamic efficiency. The 


efficiency calculated by Sage is based on the net heat 


input, including only the heat that enters the acceptor of 
the engine. These values are not directly comparable, 


and are not used for model validation, but are included 
for reference. 


Table III shows that thermodynamic power 
calculated using the current and estimated EMF is in 
good agreement with Sage predictions. 
Thermodynamic power, as measured by adding 


resistive losses to the power at the alternator leads is 


TABLE I 
MEASURED AND CALCULATED ENGINE PERFORMANCE 
Alt I*VEmMF Sage Gross Net 
Power Power Predicted Eff Eff 
+ Power 
Resistive (W) 
Loss (W) 
(W) (%) | (%) 
Meas Calc | Prediction | Meas | Sage 
Sinusoidal 435 467 457 26.2 33.2 
45 V 2nd 463 502 504 25.8 33.6 
90 V 2r 488 550 552 25.0 33.7 
110 V 2" 487 567 580 24.3 33.9 
45 V 3" 449 482 476 26.1 33.4 
90 V 3" 452 501 491 25.7 | 33.5 
135 V 3" 448 509 505 24.9 33.7 
180 V 3rd 439 530 519 23.9 33.7 
45 V 2° 
45 V 34 478 519 518 26.3 | 34.1 
90 V 2" 
45 V 3%4 484 564 573 24.7 | 33.9 
110 V 24 
45 V 34 489 579 602 23.7 | 33.9 
45 V 2" 
90 v 3% 483 541 541 25.8 | 33.9 
90 V 2" 
90 v 34 494 582 599 24.5 | 34.0 


consistently lower than the Sage predictions and the 
calculation based on EMF. The power measured at the 
terminals is considered the more direct measurement and 
is therefore used as the figure of merit in this paper. 
Figure 5 compares sinusoidal motion versus a case that 
includes the second and third harmonics. Measured 
piston position is compared to the sinusoidal case and 
instantaneous power calculated from the EMF is compared 
to the sinusoidal case and to Sage predictions. 


IV. DISCUSSION 


This experiment shows that second and third harmonic 
waveforms can be superimposed on the single harmonic 
piston waveforms typically used on Stirling engines. 
Higher harmonic waveforms with correct phasing are 
shown to increase thermodynamic power. The engines 
chosen for testing (and the vast majority of free-piston 
Stirling engine designs) have no mechanism available for 
controlling displacer motion, so optimized displacer 
waveforms could not be imposed, and displacer amplitude 
could not be limited. The resulting free-displacer motion is 
very close to sinusoidal and in many cases increases the out- 
limit of displacer motion during higher harmonic testing 
beyond the out-limit of displacer motion measured during 
the sinusoidal baseline test. Therefore, the increase in 
piston power measured during this test does not necessarily 
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Figure. 5. Comparison of sinusoidal and non-sinusoidal 
waveforms and instantaneous power. 


prove an increase in power density because the engine would have to grow to accommodate the displacer motion. 
Furthermore, since increased displacer amplitude can lead to increased power output in the absence of changes to the 
piston waveform, it is not clear how much of the power output increase measured during testing can be directly 
attributed to altering the piston waveforms and what portion should be attributed to the increase in displacer amplitude. 
However, since the displacer is a passive component driven by the pressure differential between the bounce space and 
the workings space, the increased amplitude is itself an indication of increased pressure forces and increased power 
output caused by higher harmonic piston motion. 

The thermodynamic power, calculated by several methods described in the Methodology section is in close 
agreement with analytical results provided by the Sage nodal analysis tool. However, thermodynamic power calculated 
by adding resistive alternator losses to the measured alternator power is consistently lower than calculated values and 
Sage predictions. This could be an indication of an inadequate alternator loss model or a breakdown of the assumptions 
going into the thermodynamic power calculations. Future work could address this discrepancy by developing a 
comprehensive alternator model and/or including a dynamic pressure measurement so that the thermodynamic power 
can be measured directly. In the event that future work shows Sage predictions and calculated thermodynamic power 
to be accurate, it will still be necessary to understand the loss mechanisms that prevent the thermodynamic power being 
measured at the alternator leads. 

Numerical optimization, using the Sage nodal analysis tool, identifies a second piston position harmonic of -1/2 and 
a third piston position harmonic phase angle of zero as producing the maximum increase in power output and were the 
only phase angles tested. Future work could identify other optimal waveforms based on other objective functions. 
Alternatively, a more Edisonian approach could be used in which the entire range of possible phase angles is explored, 
observing their effects on power density and efficiency. This method would eliminate the need for accurate nodal 
analysis. 

Isothermal analysis of ideal engines shows that the use of ideal piston and displacer waveforms increases Stirling 
engine power density. The magnitude of the increase depends on the method of constraint, but is found, in all cases to 
be proportional to the natural logarithm of the working space volume ratio. Isothermal analysis of a specific engine 
design including dead volume has previously shown that Stirling engines operating with ideal piston and displacer 
waveforms could increase power output 50% - 315% depending on the method of constraint as shown in [25]. Nodal 
analysis shows that loss mechanisms not taken into account in isothermal analysis reduce the power output expected 
from ideal waveforms, often times below the power output expected from sinusoidal motion. Numerical optimization 
shows that piston and displacer waveforms that compromise between the benefits of ideal motion and the reduced 
losses of sinusoidal motion can be used to increase power output by as much as 58%, depending on the method of 
constraint. Experimental data shows that engines operating with same out-limit of piston motion can be made to 
produce as much as 14% additional thermodynamic power by superimposing a second and third harmonic on the 
original sinusoidal piston waveform. 
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